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Abstract. We propose a model for the creation of entangled number states 
(Schrodinger cat states) of a Bose-Einstein condensate in a double well through simple 
phase engineering. We show that a 7r-phase imprinted condensate in a double-well 
evolves, with a simultaneous change of barrier height, to number states with well 
defined and controlled entanglement. The cat state generation is understood in terms of 
the underlying classical phase space dynamics of a 7r-phase displaced coherent state put 
at the hyperbolic fixed point of the separatrix of a physical pendulum. The extremity 
and sharpness of the final cat state is determined by the initial barrier height and the 
rate at which it is ramped during the evolution. 



The superposition principle of quantum mechanics, when applied to macroscopic 
objects, gives rise to seeming paradoxes as highlighted by Schrodinger in his famous 
cat discussion PQ. Unlike superpositions of different spin states of an electron, a 
classical object such as a room temperature cat cannot be found in a superposition 
of two macroscopically distinct states. However, carefully controlled experiments 
have been performed to create such superposition states of mesoscopic atomic and 
condensed matter systems E] • Due to the macroscopic nature of its wavefunction, 
the Bose-Einstein condensate (BEC) should be an ideal system in which to create such 
macroscopic superposition states. There have been several proposals to create cat states 
with BECs, although none have been demonstrated experimentally |U EH IE] • In 
this Letter we propose a novel model for the creation of controlled entangled number 
states of a BEC in a double well via phase imprinting on the part of the condensate 
in one of the wells followed by a change of barrier height. When properly implemented 
this results in a state of the form 



where \nL,n R ) denotes a state with n L particles in the left well, n R in the right well, 
and the total number of particles is N = + n R . Unlike in related proposals for spinor 
condensates [SHEl, we can use the barrier height to control the squeezing of the initial 
BEC ground state followed by a continuous change of barrier height to control the value 
of rii {nL = 0, 1, 2...N) of the final cat state and thereby control both the extremity (the 
value of n L ) and the sharpness (the spread around ul) of the cat state. An extreme cat 
state would correspond to rii = or N. Ground state number squeezing with variable 
barrier height in double and multi-well systems has been discussed and observed by 
many authors [HI EH HI] • 

The physical motivation for what is demonstrated here is classical: a wave-packet 
placed at the unstable fixed point of a classical phase space separatrix will bifurcate 
if allowed to time evolve. The appropriate separatrix in this case is that dividing 
bounded from rotor motion of a simple physical pendulum. It has long [T2] been 
noted that a double well quantum fluid may be modeled as a physical pendulum, an 
analysis more recently extended in Ref . fSl Ej • Modeling a repulsive BEC in a double 
well in the two mode approximation jHJ E3, the coefficients of the exact Fock space 
eigenstates may be projected into a semi-classical (n, 9) phase space via the Husimi 




(1) 
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probability distribution [ElUZj, whereby the pendulum analogy is immediately seen to 
be a property of the exact quantum eigenfunctions, as illustrated in Fig. The low 
energy states correspond to oscillatory motion of the pendulum, and the higher-lying 
states to superpositions of pendulum rotor motions in two opposite directions, these 
being macroscopic superposition states. The ground state in this model is a Gaussian 
centered at the origin in (n, 6) space with number squeezing controlled by the barrier 
height 0. A displacement in phase by n, as may be created by phase shifting the single 
particle wavefunction in one of the wells, shifts the ground state wave-packet to the 
unstable fixed point. As expected, the subsequent autonomous dynamics of the shifted 
packet lead to bifurcation and production of macroscopic superposition states. These 
may be further controlled, leading to sharp number entangled states with a chosen value 
of til, as in Eqn. (1), by choice of the initial state and barrier ramping, as discussed 
below. Further details of the quantum phase space picture and its applications are 
presented elsewhere [TH] . 

The many-body Hamiltonian for a system of iV weakly interacting bosons in an 
external potential V(r) in second quantization, is given by 

h 2 

•7^V 2 + \/(r) 



H = / dr&( 



2m 

9 



+| J dr&(r)&(r)$(r)$ (r) (2) 

where ^(r) and ^(r) are the bosonic annihilation and creation field operators, m is the 
particle mass, and g = Anc ^ h where a s is the s-wave scattering length. 

Many authors have studied the double-well condensate using the two-mode 
approximation P [T3J . We use the model introduced by Spekkens and Sipe 0. The 
exclusion of the nonlinear tunneling terms in this model gives rise to the Bose-Hubbard 
model JH] • The full two-mode Hamiltonian is 

H = e LL N L + srrNr + (e LR + gT^N - l))(a{a R + a R a L ) 
9T2 

— ^(aL a L a R a R + a R a R a L a L + 4:N L N R ) (3) 



where Nl 




(4) 



where i,j=L,R. 

T = J dr^Kr);^ = J dr0|(r)^(r); T 2 = J dr^(r)4(r). (5) 

€ll and e RR are the energies of a single particle in the left and right wells, each is 
multiplied by the number operators Nl,r counting the number of particles in each well; 
€l R is the single particle tunneling amplitude which appears as a coefficient of operators 
that allow a single particle to hop from one well to the other; T is the nonlinear mean- 
field energy in each well and T\ 2 are nonlinear tunneling matrix elements. 
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To simplify a theoretical study of this two-mode Hamiltonian, we make a one 
parameter approximation of the single particle energies and the tunneling matrix 
elements: 

£ll = e RR = T = 1; e LR = T x = -e' a ; T 2 = -e~ 2 °. (6) 

This parametrization allows a simple study of continuous change in the linear and non- 
linear tunneling through variation of a single parameter a, which is associated with the 
barrier height and width. In our computations with this model we ignore the T 2 term 
which scales as e~ 2a . 

Taking into account the conservation of the total number of particles, the state space 
for the system can be spanned in iV + 1 Fock state basis \ni, N — ul), til = 0, 1, ...N, 
where ul particles occupy the single particle state in the left well and N — til particles 
occupy the right well. The state vector is a superposition of all the number states 

N 

l*>= E 4\\n L ,N-n L ) (7) 



n L =0 



(J \n L (J \N~n L 

\n L , N-n L )= {aL l_ {a £ \vac). (8) 



where 



Finding the eigenvalues and eigenvectors of the model Hamiltonian in the Fock basis 
can be easily accomplished by diagonalizing a (N + 1) x (TV + 1) tridiagonal matrix. 
Fig. |21 panels (a) and (b) show the coefficients of the eigenvectors for the two lowest 
lying states for 40 particles for a = 4. Panels (c) and (d) show the 30th and the 31st 
eigenstates. The lowest lying states clearly have harmonic oscillator like wavefunctions 
in Fock space. On the other hand the higher lying states are cat-like, i.e. superpositions 
of states with most particles in the left and right wells. 

The characteristics of these states can be understood in terms of the motion of 
a momentum-shortened physical pendulum. Anderson suggested that the dynamics 
of the Josephson effect can be understood in terms of a physical pendulum [T2] ; 
similarly the mean-field dynamics of a Bose-Einstein condensate in a double well is well 
described by a momentum-shortened physical pendulum, which has novel phase space 
characteristics [TBI EH EH ED) • It is natural to ask what aspects of the mean-field phase 
space properties are contained in the exact quantum treatment. Phase space formulation 
of quantum mechanics is extremely useful in studying many aspects of quantum optics, 
collision theory and quantum chaos, and this is the approach we take to study the 
double well. The Husimi probability distribution [TBI El can be used to project, in 
a squeezed coherent state representation, the classical phase space properties from 
stationary quantum wavefunctions. For a BEC in a double well, the phase difference 
6=6l-6r and the number difference n = nL ~ nR are the conjugate variables analogous to 
q and p. In this (n,0) representation a probability distribution function can be defined 
as 

P J (n,6) = \(9 + in\y i )\ 2 (9) 
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where 

1 N/2 , , _ n2 

^ + m l^ = 7^ E 4'exp[^n'-^— % (10) 

Here n' = nL ~ nR , rather than being the simpler left particle counter, and c n > are the 
corresponding Fock-state coefficients. The Q-function in quantum optics is a special case 
of the Husimi distribution function where k=u, u being the frequency of a coherent state 
Gaussian wave-packet ^Zj. The 'coarse- graining' parameter k determines the relative 
resolution in the conjugate number-phase (n,#) space. 

Fig. [T] panels (b), (c), (d) and (e) show the density plots for Husimi probability 
distributions for 40 particles for the ground state, 6th, 12th and 35th states respectively. 
Panel (a) shows the classical energy contours that are also the classical trajectories. As 
is evident from the panels, the ground state is a minimum uncertainty wave-packet in 
both number and phase that is centered at the origin, the harmonic-oscillator-like low 
lying excited states are the analog of pendulum librations, and the higher lying cat-like 
states are the analog of rotor motions, with a clear signature of the quantum separatrix 
state where the libration and rotation states separate. Note that the Husimi distribution 
of the 35th state shown in panel (e) is a superposition of most particles in the left and 
right wells. In the classical phase space this corresponds to two trajectories above and 
below the separatrix that denote rotor motions of a physical pendulum in two opposite 
directions. The creation of such macroscopic superposition states with BEC is possible 
with simple phase engineering. It is clear that near classical behavior is observed for 
quite small values of N. Larger N further sharpens the classical correspondence. 

Writing phases on part of a condensate is experimentally feasible via interaction 
with a far off-resonance laser. This method has been used to generate dark solitons 
and measure their velocities due to a phase differential |2~Tj |2"2"] . Mathematically, such a 
method corresponds to multiplying the coefficient of each of the Fock state basis of an 
eigenstate with e~ mLe , where \til) is the corresponding Fock state, and 6 is the phase 
offset for particles in the left well. By 7r phase imprinting the condensate in one well, 
the ground state centered at the origin (0,0) in phase space is displaced to the unstable 
equilibrium point (0,7r) on the separatrix. Using exact quantum time evolution within 
the framework of the two mode model, the resulting quantum wave-packet bifurcates 
as expected. If the barrier is raised as discussed below, the wave-packet is permanantly 
split, resulting in a superposition of two classical rotor states. With the appropriate 
ramping of the potential barrier and proper tuning of the values of the parameters in 
the model we can vary the sharpness and extremity of the cat states. 

As an example, in Figs. El and 0] we show how a cat state with 1000 particles is 
generated. Fig. El shows the evolution in phase space using Husimi projections - (a) the 
ground state, (b) a 7r-phase imprinted state, (c) and (d) show subsequent evolution in 
the process of bifurcating the state; further evolution along with a change of barrier 
totally splits and traps the state above and below the separatrix, as shown in (e), 
finally giving rise to cat state in (f). Shown in Fig. H]is the evolution in the Fock state 
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basis. A close view of the coefficients, shown in (f), shows that it is a rather sharply 
peaked entangled number state. The parameters used for this example are g — Tq — 1, 
€lr = T\ = —e~ a where the barrier height is ramped up in time as a = 3 + 2t. When a 
cat state is reached the barrier is suddenly raised to essentially halt the evolution. With 
different initial barrier heights and the same ramping of the potential, the extremity of 
the cat states can be tuned. Examples are shown in Fig. |3] where the different values 
of the barrier heights are a = 1 + 2t, a = 3 + 2t and a = 5 + 2t for rows (1), (2) and 
(3) respectively. The columns show: (a) the barrier height and the ramping, (b) the 
respective ground state, (c) the final cat state at the end of the ramping, and (d) a close 
view of the coefficients for the final state. As is evident from the pictures, the initial 
squeezing of the ground state determines the extremity of the final cat state. The rate 
at which the barrier is ramped determines the sharpness. It is important to be able to 
tune to less extreme cat states as such states are more robust to loss and decoherence. 

In conclusion, we have shown that a 7r-phase imprinted BEC in a double well evolves 
to an entangled number state (Schrodinger cat state) which is a superposition of most 
particles being in the left and right wells simultaneously. Initial squeezing of the ground 
state and a ramping of the barrier height help tune the extremity and sharpness of the 
final cat state. A semi-classical phase space analysis of the full quantum problem reveals 
the similarity with a classical pendulum phase space, thus providing a visual explanation 
of cat state generation through motions along the separatrix. Making condensates in 
a double-well, phase imprinting and tuning the barrier are viable current experimental 
techniques. However, decoherence mechanisms will put severe constraints on the size or 
extremity of cat states that can be produced in a laboratory. Such challenges will make 
the attainment of cat states with a BEC interesting and important. 

We would like to thank Sarah B. McKinney for discussions and computational 
support. This work was supported by NSF grant PHY-0I4009I. 
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Figure 1. Comparison of the classical nonrigid pendulum phase space and Husimi 
probability distributions for different energy eigenstates for 40 particles. Shown arc 
(a) classical energy contours. Husimi projections for (b) the ground state which is 
a minimum uncertainty wave-packet centered at the origin, (c) the 6th state which 
is harmonic-oscillator like and the analog of pendulum librations, (d) the 12th state 
which is a quantum separatrix state separating the libration and rotor states, (e) the 
35th state which is analogous to a superposition of classical pendulum rotor motions 
in two opposite directions. The creation of such macroscopic superposition states with 
a BEC is possible with phase engineering. 
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Figure 3. Shown is the evolution to a cat state in Husimi projection space, (a) The 
ground state at t=0, (b) the 7r-phase imprinted ground state at the hyperbolic fixed 
point, (c) at t=0.01 the wave-packet is bifurcating along the separatrix, (d) at t=0.016 
it continues to move along the separatrix, (e) at t=0.4 the states become trapped as 
we increase the barrier, and (f) at t=2.3 a sharply peaked cat state is obtained. 
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Figure 4. Shown is the evolution to a cat state in Fock space at the same time 
instants as the previous figure: (a) t=0, (b) t=0.01, (c) t=0.016, (d) t=0.4, (e) t=2.3, 
(f) is a magnified version of (e), showing the nonvanishing Fock state coefficients. Note 
that for clarity the probability amplitudes in the vertical axes have different scalings 
for the panels. 
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Figure 5. Shown are the cat states with different initial heights of the barrier and 
therefore different initial squeezings of the BEC ground state, but the same ramping 
of the potential. Row (1) shows the states where a = 1 + 2t: (a) the parameter a as 
a function of time, (b) the ground state, (c) the final cat state, and (d) a magnified 
view of the Fock-state coefficients. Rows (2) and (3) show the results for a = 3 + It 
and a — 5 + 2t respectively. The initial barrier height controls the extremity of the 
cat states. Note that for clarity the axes in the panels have different scalings. 



